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Abstract 

The recent discoveries of novel human coronaviruses, including the coronavirus causing 
SARS, and the previously unrecognized human coronaviruses HCoV-NL63 and HCoV- 
HKU1, indicate that the family Coronaviridae harbors more members than was previously 
assumed. All human coronaviruses characterized at present are associated with respira¬ 
tory illnesses, ranging from mild common colds to more severe lower respiratory tract 
infections. Since the etiology of a relatively large percentage of respiratory tract diseases 
remains unidentified, it is possible that for a certain number of these illnesses, a yet 
unknown viral causative agent may be found. Screening for the presence of novel coro¬ 
naviruses requires the use of a method that can detect all coronaviruses known at present. 
In this chapter, we describe a pancoronavirus degenerate primer-based method that allows 
the detection of all known and possibly unknown coronaviruses by RT-PCR amplification 
and sequencing of a 251-bp fragment of the coronavirus polymerase gene. 

Key words: coronavirus; pancoronavirus RT-PCR; degenerate primers; sequencing; 
polymerase gene 


1. Introduction 

At present, viral culture is the “gold standard” for laboratory diagnosis of 
respiratory infections. Since coronaviruses are very difficult to grow in cell 
culture, accurate and sensitive diagnoses are not feasible by this technique. To 
overcome the lack of sensitivity and to obtain rapid diagnostic results, more 
sensitive molecular methods for the detection of human coronaviruses (HCoVs) 
have been developed, including reverse-transcriptase polymerase chain reaction 


From: Methods in Molecular Biology, vol. 454: SARS- and Other Coronaviruses , 

Edited by: D. Cavanagh, DOI: 10.1 007/978-1-59745-181-9_1, © Humana Press, New York, NY 


3 



4 


Vijgen et al. 


(RT-PCR), nested RT-PCR, and recently real-time RT-PCR (7). Nevertheless, 
HCoVs are not often diagnosed in clinical laboratories, although the SARS epi¬ 
demic and the identification of HCoV-NL63 drew attention to the clinical rele¬ 
vance of HCoVs. Multiplex RT-PCRs for the detection of common respiratory 
viruses in clinical specimens have now been supplemented with primer pairs for 
amplification of HCoV-OC43, 229E, and NL63 (2,3). 

In order to allow the detection of all known coronaviruses in one assay, 
consensus RT-PCRs have been developed, based on an alignment of conserved 
genome regions of several coronaviruses (4). However, the recent identification 
of HCoV-NL63 demonstrated that this novel coronavirus could not be amplified 
by the consensus RT-PCR described by Stephensen and colleagues, which might 
be explained by the presence of several mismatches in the primer sequences. We 
designed a novel pancoronavirus RT-PCR in which we modified the coronavirus 
consensus RT-PCR primers based on an alignment of the HCoV-NL63 sequence 
and the corresponding sequences of 13 other coronaviruses. Theoretically, this 
pancoronavirus RT-PCR should amplify a 251-bp fragment of the polymerase 
gene of all coronaviruses known at present, and we tested this experimentally 
for the five known human coronaviruses (HCoV-NL63, 229E, OC43, HKU1, 
and SARS-CoV) and three animal coronaviruses (PHEV, FIPV, and MHV). The 
sensitivity of the pancoronavirus RT-PCR was found to be lower than for a non¬ 
degenerate RT-PCR, which can be explained by the high level of degeneracy of 
the pancoronavirus primers. Nevertheless, this pancoronavirus assay is a useful 
tool for screening sample collections for the presence of all known and poten¬ 
tially yet unknown coronaviruses (5). 


2. Materials 

2.1. Sample Handling and RNA Extraction 

1. Viral transport medium (VTM) (BD Diagnostic systems, USA) or phosphate 
buffered solution (PBS). Store at room temperature. 

2. Filter Minisart plus 0.45 pm (Sartorius, Vivascience AG). 

3. 2-ml syringe without a needle. 

4. Glass pearls (5 mm diameter) (Merck MDA, VWR, Belgium). 

5. Qiamp viral RNA mini kit (Qiagen, The Netherlands). Store at room temperature, 
except for the AVL lysis buffer, which should be stored at 4°C after addition of 
carrier RNA. 


2.2. Pancoronavirus One-Step RT-PCR 

1. Primers: Forward primer: Cor-FW (5 / -ACWCARHTVAAYYTNAARTAYGC-3 / ) 
and reverse primer: Cor-RV (5 / -TCRCAYTTDGGRTARTCCCA-3 / ) (Eurogentec, 
Seraing, Belgium). The primers are prepared in a stock solution of 100 pM. 
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A working solution of 15 p,M is made for use in the RT-PCR. Store at — 20°C 
until use. 

2. Qiagen One-Step RT-PCR kit containing One-Step RT-PCR enzyme mix (a com¬ 
bination of Omniscript and Sensiscript reverse transcriptase and HotStarTaq DNA 
polymerase), One-Step RT-PCR buffer 5X, and dNTP mix, 10 mM each (Qiagen, 
The Netherlands). Store all reagents at — 20°C. 

3. RNAse-free water (Sigma-Aldrich, Belgium). Store at room temperature. 

4. GeneAmp PCR system 9700 thermal cycler (Applied Biosystems, Foster City, 
CA, USA). 

5. IX Tris-borate EDTA (TBE) buffer: dilute 10X TBE buffer (Invitrogen Life Tech¬ 
nologies) 1:10 with MilliQ water (Millipore). Store at room temperature. 

6. PAGE gels (6% acrylamide/N^-methylenebisacrylamide in IX TBE buffer con¬ 
taining 0.064% (w/v) ammonium persulfate and 150 pd tetramethylethylenedi- 
amine). After preparation of a batch of ten gels store at 4°C. 

7. Molecular Weight Marker VI (MWM VI) (Roche, Mannheim, Germany). Store 
at 4°C. 

8. Imagemaster VDS (Pharmacia, Uppsala, Sweden). 


2.3. Sequencing and Sequence Analysis 

1. MSB Spin PCRapace purification kit (Invitek, Westburg, The Netherlands). Store 
at room temperature. 

2. Saekem GTG agarose for the recovery of nucleic acids (Cambrex, Belgium). 

3. IX Tris acetate EDTA (TAE) buffer: dilute 10X TAE buffer 1:10 with MilliQ 
water (Millipore). Store at room temperature. 

4. QIAquick gel extraction kit (Qiagen, The Netherlands). 

5. Forward primer: Cor-FW (5 / -ACWCARHTVAAYYTNAARTAYGC-3 / ) or 
reverse primer Cor-RV (5 / -TCRCAYTTDGGRTARTCCCA-3 / ) (Eurogentec, 
Seraing, Belgium). The primers are prepared in a stock solution of 100 pM. A 
working solution of 5 pM is made for use in the cycle sequencing reaction. Store 
at — 20°C until use. 

6. ABI PRISM BigDye Termination Cycle Sequencing Ready Reaction kit (Applied 
Biosystems, CA, USA). Store at — 20°C. 

7. Biometra T3000 thermocycler (Biometra, Westburg, The Netherlands). 

8. Ethanol absolute. 

9. Sodium acetate 3 M pH 4.6. Store at 4°C. 

10. Formamide. Store at — 20°C. 

11. ABI Prism 3100 Genetic Analyzer (Applied Biosystems, Foster City, CA, USA). 

3. Methods 

3.1. Sample Handling and RNA Extraction 

1. Nasopharyngeal and throat swabs: place the tip of the swab in a tube containing 
3 ml VTM or PBS. After squeezing the tip of the swap and vortexing, take 140 pd 
of the solution for RNA extraction. Store the rest of the sample at — 80°C. 
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2. Bronchoalveolar lavages and nasopharyngeal aspirates: take 140 pd of the sample 
for RNA extraction without prior dilution in VTM or PBS. Store the rest of the 
sample at — 80° C. 

3. Sputum samples: dilute about 100 to 300 mg of sputum in 3 ml VTM or PBS. 
Add glass pearls and vortex intensively until the solution containing the sputum is 
completely liquefied. Filter the solution using a 2-ml syringe without a needle and 
a 0.45-jam Minisart Plus filter (see Note 1). Take 140 pd of the filtrate for RNA 
extraction. Store the rest of the sample at — 80°C. 

4. Perform the RNA extraction on 140 pd of sample according to the manufacturer’s 
instructions using the Spin Protocol (Qiagen QIAamp viral RNA mini kit hand¬ 
book) (see Note 2). Store the obtained RNA extract (60 pd) at — 80°C. 


3.2. Pancoronavirus One-Step RT-PCR 

1. Screening a sample for the presence of coronaviruses is performed by amplify¬ 
ing a 251-bp fragment of the coronavirus polymerase gene using the following 
primer set: Cor-FW (5 / -ACWCARHTVAAYYTNAARTAYGC-3 / ) and Cor-RV 
^-TCRCAYTTDGGRTA RTCCCA-30 (see Fig. 1) (see Note 3). Prepare the 
following master mix for RT-PCR according to the number of samples to be tested: 
10 [jlI 5X QIAGEN OneStep RT-PCR Buffer, 2 pd dNTP mix (final concentration 
of 400 pM of each dNTP), 1.8 pd QIAGEN OneStep RT-PCR Enzyme Mix, 4 p,M 
of forward and reverse primer, and RNase-free water to 40 pd per reaction (see 
Notes 4 and 5). 


HCOV-NL63 
UCoV-2 2 9E 
FI PV 
TGSV 
PFDV 

HCoV-OC43 

BCoV 

PHSV 

CRCV 

MHV 
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ACWCARHTVAAYYTNAARTAYGC 


TGGGAYTAYCCHAARTGYGA 


Fig. 1. Selection of primers for the novel pancoronavirus RT-PCR. Shown is the 
alignment of 14 coronaviral sequences of a conserved region of the polymerase gene. 
The forward (Cor-FW) and reverse (Cor-RV) primer sequences are shown at the bottom 
(Y=C/T, W=A/T, V=A/C/G, R=A/G, H=A/T/C, N=A/C/T/G). 
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2. Add 10 fjul RNA-extract to the tube containing 40 jjlI of the master mix, leading 
to a final reaction volume of 50 jjlI (see Notes 6 and 7). The reaction is carried 
out with an initial reverse transcription step at 50°C for 30 min, followed by PCR 
activation at 95°C for 15 min, 50 cycles of amplification (30 sec at 94°C; 30 sec at 
48°C; 1 min at 72°C), and a final extension step at 72°C for 10 min in a GeneAmp 
PCR system 9700 thermal cycler (Applied Biosystems, Foster City, CA, USA). 

3. Load 9 jxl of PCR product mixed with 1 pi loading dye on a 6% PAGE gel. Load 
3 pi of MWM VI on each side of the gel (see Note 8). Run the gel at 200 V, 
200 mA for 30 min. Visualize the bands by illuminating the gel with UV-light 
on an Imagemaster VDS system after staining the gel for 5 min with ethidium 
bromide (EtBr) (see Note 9). 

4. Determination of the fragment size is performed by comparing the length of the 
bands to MWM VI. Samples from which a PCR product of approximately 251 bp 
is amplified are presumed to be coronavirus-positive and will be subject to further 
(sequence) analysis (see Note 10). 


3.3. Sequencing and Sequence Analysis 

1. If a single 251-bp band is visible after EtBr staining, purify the total ampli¬ 
fication product by using the MSB Spin PCRapace purification kit according 
to the manufacturer’s instructions (Invitek MSB Spin PCRapace purification kit 
handbook). 

2. If more than one band is visible, including a 251-bp band, run the total amount 
of PCR product mixed with 5 pi loading dye on a 1.5% agarose gel (in IX TAE 
buffer) on 100 V, 500 mA for 1 h. Stain for 10 min with EtBr, and cut each band 
from the gel with a sterile scalpel. Purify the fragments excised from the agarose 
gel using the QIAquick gel extraction kit according to the Spin Protocol in the 
QIAquick gel extraction kit handbook. 

3. For sequencing, add 5 pi of purified amplification product to 4 pi of ABI PRISM 
BigDye Ready Reaction mix and 1 pi of a 5-pM solution of forward or reverse 
primer. The sequence reaction, which consists of 25 cycles of 30 sec at 96°C, 15 
sec at 50°C, and 4 min at 60°C, is performed in a Biometra T3000 Thermocycler. 

4. After cycle sequencing, precipitate the DNA products by adding 62.5 pi absolute 
ethanol, 3 pi sodium acetate, and 24.5 pi MilliQ water per sample. Mix the sam¬ 
ples and leave them at room temperature for 15 min, followed by a centrifugation 
for 20 min at 13,000 rpm (see Note 11). Carefully aspirate the supernatant, wash 
the pellet by adding 150 pi 70% ethanol, and centrifuge the samples for 5 min at 
13,000 rpm. Carefully aspirate the supernatant, and let the pellet dry for 15 min 
at room temperature. Dissolve the pellet in 15 pi of formamide and denature for 
2 min at 95°C. The samples are then stored at —20°C before analysis on an ABI 
Prism 3100 Genetic Analyzer (Applied Biosystems, Foster City, CA, USA). 

5. Chromatrogram sequencing files can be inspected with Chromas 2.2 (Technely- 
sium, Helens vale, Australia), and contigs of the fragments sequenced by Cor-FW 
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and Cor-RV, respectively, can be prepared using SeqMan II (DNASTAR, Madison, 
WI). The contig sequence data can easily be arranged into a FASTA format by 
saving the consensus as a single text file in the Seqman program. The text file will 
automatically have the following format: 

> sequence title (hard return) 

GACGAGTAATTGC... (not containing hard returns) 

Such a FASTA file can then be used for nucleotide similarity searches using 
the NCBIWWW-BLAST (basic local alignment search tool) server on the Gen- 
Bank DNA database release 155.0 ( 6 ). The coronavirus type that is present in the 
sample can be determined by using this BLAST analysis, which provides coron- 
avirus sequences that show the highest nucleotide similarity with the sequenced 
fragment (see Note 12). Furthermore, a FASTA file of the sequence data can 
be loaded into a multiple sequence alignment program such as CLUSTAL X 
version 1.82 (7), and the resulting alignment can be manually edited in the 
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Fig. 2. Neighbor-joining phylogenetic tree of coronavirus partial polymerase gene 
sequence data (corresponding to the pancoronavirus fragment). The frequency of occur¬ 
rence of particular bifurcations (percentage of 500 bootstrap replicate calculations) 
is indicated at the nodes. Bootstrap values over 75% are shown. The coronavirus 
sequences used here are available from GenBank under the following accession num¬ 
bers: HCoV-NL63, AY567487; HCoV-229E, AF304460; PEDV, AF353511; transmissi¬ 
ble gastroenteritis virus (TGEV), AF304460; feline infectious peritonitis virus (FIPV), 
AF124987; HCoV-OC43, AY391777; HCoV-HKUl, NC 006577; porcine hemagglu- 
tinating encephalomyelitis virus (PHEV), AF 124988; bovine coronavirus (BCoV), 
AF391541; mouse hepatitis virus (MHV-A59), X51939; sialodacryoadenitis virus 
(SDAV), AF124990; canine respiratory coronavirus (CRCV), AY150273; SARS-CoV, 
AY313906; turkey coronavirus (TCoV), AF124991; and infectious bronchitis virus 
(IBV), Z30541. 
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GeneDoc Alignment editor (8). Starting from a sequence alignment that com¬ 
pares nucleotides at homologous sites, i.e., sharing a common ancestor, a phylo¬ 
genetic tree can be constructed using phylogenetic programs such as MEGA ver¬ 
sion 3.1 (Fig. 2) (9). The construction of a phylogenetic tree based on sequence 
data of samples and coronavirus sequence data from GenBank will provide 
information regarding the genetic and evolutionary relatedness of a coronavirus 
detected in a sample to other known coronaviruses (see Note 13). 


4. Notes 

1. Sputum samples diluted in VTM or PBS should be filtered prior to analysis in 
order to remove large amounts of cellular debris. It can be difficult to completely 
liquefy the sputum sample, so the filter may become blocked. In this case, the 
only solution is to use another filter. 

2. Especially during RNA extraction, the likelihood of cross-contamination is high. 
When performing the RNA extraction manually, one should change tips in 
between samples, as well as when adding wash buffer to the samples. Another 
precautionary measure that can be taken is to use tissue paper when opening the 
samples to avoid the formation of aerosols. 

3. For modification of the coronavirus consensus RT-PCR primers (4), an align¬ 
ment of the polymerase gene sequences of 14 known coronaviruses was made 
(Fig. 1). Degenerate Cor-FW and Cor-RV primers were developed. The coordi¬ 
nates of Cor-FW and Cor-RV are 14017 and 14248, respectively, in the HCoV- 
NL63 complete genome sequence. The 14 coronavirus sequences used here 
are available from GenBank under the following accession numbers: HCoV- 
NL63, AY567487; HCoV-229E, AF304460; PEDV, AF353511; transmissible 
gastroenteritis virus (TGEV), AF304460; feline infectious peritonitis virus 
(FIPV), AF124987; HCoV-OC43, AY391777; porcine hemagglutinating 
encephalomyelitis virus (PHEV), AF124988; bovine coronavirus (BCoV), 
AF391541; mouse hepatitis virus (MHV-A59), X51939; sialodacryoadenitis 
virus (SDAV), AF124990; canine respiratory coronavirus (CRCV), AY150273; 
SARS-CoV, AY313906; turkey coronavirus (TCoV), AF124991; and infectious 
bronchitis virus (IBV), Z30541. 

4. The sensitivity of the pancoronavirus RT-PCR assay was assessed by testing 
tenfold dilutions of HCoV-NL63 and HCoV-OC43 RNA. Primer concentrations 
were increased up to 4 pM per reaction to improve the sensitivity of the assay. 
While 50 copies of HCoV-OC43 RNA copies per microliter sample could be 
detected, the sensitivity for HCoV-NL63 was slightly lower, i.e., 5 x 10 3 RNA 
copies per microliter sample. 

5. An important precautionary measure when performing RT-PCR is to avoid 
RNase contamination. Therefore, all materials should be RNase-free, and gloves 
and a gown should be worn even when preparing the (noninfectious) master mix. 
In addition, an RNase inhibitor (RNase OUT, Invitrogen Life Technologies; 20 
units per reaction) can be added to the RT-PCR reaction to inhibit degradation 
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of the RNA sample owing to RNases. This measure should be taken especially 
when a positive control RNA sample no longer amplifies. 

6. Include negative as well as positive controls in each RT-PCR. A negative control 
consists of a sample containing only the master mix without RNA (10 jxl of water 
is added instead of RNA). As a positive control, RNA of a human coronavirus 
(e.g., HCoV-OC43) can be used. 

7. When adding the RNA samples to the master mix, precaution must be taken to 
avoid RNA contamination. Such a precautionary measure is opening the tubes 
with tissue paper, thereby avoiding the formation of aerosols, as the latter could 
lead to cross-sample contamination. 

8. It is advised to load the Molecular Weight Marker on both sides of the gel, as a 
gel can show the phenomenon of “smiling,” which makes the determination of 
the fragment size difficult when the MWM is only present on one side of the gel. 

9. Use two pairs of gloves when handling ethidium bromide, as this intercalating 
agent is highly mutagenic. 

10. Amplification of a 251-bp fragment of the coronavirus polymerase gene was 
demonstrated for HCoV-NL63, HCoV-OC43, HCoV-229E, SARS-CoV, PHEV, 
FIPV, and MHV-A59 (Fig. 3) and for HCoV-HKUl (data not shown). 

11. The DNA pellet is often not visible. Therefore, during centrifugation the tubes 
are always put in the same orientation in the centrifuge, so that one knows at 
which side of the tube the pellet is located. The aspiration of supernatant should 
be performed carefully without touching the pellet. 

12. The polymerase gene is a very conserved region in the coronavirus genome. 
Therefore, the percentage sequence similarity will be high among different coro¬ 
navirus members belonging to the same group. Nevertheless, this sequence infor¬ 
mation allows a primary identification of the coronavirus type that is present in 
a sample. 



Fig. 3. Gel electrophoresis after pancoronavirus RT-PCR assay. The indicated band 
of 251 bp corresponds with the expected amplicon size. Molecular Weight Marker VI 
was used as a marker (Roche, Mannheim, Germany). 
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13. The genetic and evolutionary relatedness to known coronaviruses can be deter¬ 
mined from the phylogenetic tree clustering pattern of the coronavirus sequence 
fragment, amplified from a sample. A phylogenetic tree is composed of nodes 
that are connected by branches. The external nodes represent the existing organ¬ 
isms or taxa, while the internal nodes represent the (hypothetical) ancestors of 
the existing taxa. The length of the horizontal branches is indicative of the evo¬ 
lutionary distance between the connected nodes, whereas the vertical branches 
are noninformative. Taxa that are present in the same cluster are monophyletic, 
which indicates that they show a (relatively) close genetic relatedness. 
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